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Point field models for the galaxy point pattern 
Modelling the singularity of the two-point correlation function 
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Abstract. There is empirical evidence that the two-point correlation function of the galaxy distribution follows, 
for small scales, reasonably well a power-law expression f(r) oc r~ 7 with 7 between 1.5 and 1.9. Nevertheless, 
most of the point field models suggested in the literature do not have this property. 

This paper presents a new class of models, which is produced by modifying point fields commonly used in 
cosmology to mimic the galaxy distribution, but where 7 = 2 is too large. The points are independently and 
randomly shifted, leading to the desired reduction of the value of 7. 
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1. Introduction 



Saar 20*02] ). It is also worth mentioning that two point 



The two-point correlation function £(r) has been the pri- 
mary tool for quantifying large-scale cosmic structure (see 
Peebles 1980 ). Several estimates of this statistical quantity 
from rcdshift surveys suggest a power-law shape of £(r) for 
0.5h -1 Mpc < r < 10h _1 Mpc (h being the Hubble con- 
stant in units of 100 km s" 1 Mpc -1 ) 



£(r) 



(1) 



with 7 betwee n 1.5 and 1.9, sec for e xample |Davis fc 



Peebles (1983 


)■ 


Maddox et al. (199C|), 


Martinez et al. 


(1993|) or Pons-Borden'a et al. (1999|). Thus, £(r) seems 



to have a singularity of order 7. However, taking into ac- 
count that galaxies have a physically finite extended size, 
the shape of £(r) has necessarily to deviate, at very small 
scales, from the power-law behavior. 

It is important to be able to generate controllable 
point processes resembling the galaxy distribution, not 
only showing visual similar patterns to the observed one, 
but also reproducing the same statistical properties. These 
processes are necessary to test the performance and, in 
particular, the discriminative power of the statistical es- 



fields could present similar second-order characteristics, 
such as the two-point correlation function analyzed in this 
paper, but show completely different results for higher or- 
der statistics. 

Several methods have been published so far to gener- 
ate mock catalogues to mimic the galaxy distribution. We 
could distinguish two main kinds of methods: Those based 
on a physically motivated procedure and those based on a 
stochastic model with controllable statistical properties. 
The first family of methods includes the following two 
main approaches: 

1. Starting with cosmological A-body simulations based 
on the gas-dark matter dynamics, one can extract 
galaxy catalogues by applying physical schemes related 



with the galaxy formation processes (Kauffmann et al 



timators used to describe the cosmic texture (Martinez 



1999 ). This is obviously the most realistic way to ob- 
tain mock representations of the galaxy distribution. 
However, it requires much computing time and there- 
fore the cost of obtaining several realizations is high. 
The second possibility is to treat the dynamics of 
the non-linear clustering of matter by means of ap- 
proximate methods based on higher order perturba- 
tion theory. A recent example is the so-called PThalo 



model introduced by Scoccimarro & Shcth (2002) 
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These authors select virialized dark matter halos from 
simulations obtained by applying the second-order 
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Lagrangian perturbation theory. Then, given a bias- 
ing scheme which accounts for the fact that the galaxy 
formation needs not to be proportional to the local 
matter density, one can extract realistic galaxy distri- 
butions. 

The second family of models for the galaxy distribution 
are those which are based on statistical or geometrical pre- 
scriptions. These models are more flexible than those from 
the first family, can be easily simulated and manipulated, 
and have more controllable statistical properties. In many 
cases, these statistics can be calculated analytically. The 



first model of this family was introduced by Neyman & 



Scott ( 195S ) and became very popular in spatial statistics. 



The characteristics of this model w ere probably prefigured 
by the Lick survey. Peebles (198d| ) reviews two point pro- 
cesses providing clustering patterns with power-law be- 
havior for the two-point correlati on function: The Soneira 
and Peebles hierarchical model ( Boneira fc Peebles 1978 ) 
and the Mandelbrot' s prescription base d on the Rayleigh- 
Levy random walk ( Mandelbrot 1982 ). The first model, 
based on a superposition of several fractal clumps, repro- 
duces well the second-, third-, and fourth-order correlation 
functions. In the Rayleigh-Levy model, galaxies are placed 
at the steps of a random walk starting from an initial po- 
sition, where each jump takes an isotropically chosen ran- 
dom direction with the length following a power-law prob- 
ability distribution function. Although this model shows 
point distributions with little resemblance to the observed 
one, it has been successfully applied to test the error mea- 



surements of counts- in-cells statistics ( Szapudi & Colombi 
1996|)~|The singularity of the correlation function in these 
models can be varied in a rather wide interval. 

Other point field models published in literature present 
a two-point correlation function without any singularity 
or with a singularity of order 2. A well-known example 
is the Cox segment process, where points are situated on 
randomly orien ted line segments. I ts correlation function 
was derived by Stoyan et al. (1995 ) and has a singularity 
of order 2. A similar cluster process was introduced by 
Buryak fc Doroshkevich (1996| ), where the points (i. e., the 
galaxies) populate a mixture of randomly oriented lines 
and planes. The singularity of the correlation function is 
determined by the line component and has again the order 
2. 

Van de Weygaert (1991) suggests a point field model 
which is built from the vertices of a Voronoi tessella- 
tion calculated from a randomly distributed nuclei. This 
model, based on one of the most interesting constructions 
studied in stochastic geometry, has the property of mim- 
icking the statistical properties of the distribution of clus- 
ters of galaxies. However, for this model, the order of the 
singularity of is again 2. 



A still simpler point field model is appearing in 3toyan 



(1992). It is a Poisson cluster process where each cluster 



consists of only two points. The corresponding two-point 
correlation function can be given a singularity of any de- 
sired order between and 3. Of course, this point field is 



rather primitive and has no physical motivation. It makes 
sense only with regard to the singularity of the two-point 
correlation function. 

This paper investigates an idea which leads to a new 
class of point field models. They can have any desired 
value of 7 between and 3 and seem to be more realistic 
than the previous models. 

The models listed above are based on simplifying as- 
sumptions which are not well fulfilled in reality. The idea is 
now to model the difference between ideal and real world 
by means of an additional random element, i. e., by giving 
independently a random shift to each point of the theo- 
retical model. 

The analytical investigation of the effect of these ran- 
dom shifts is very difficult and is therefore carried out only 
for the case of relatively simple Poisson cluster processes, 
the Neyman-Scott processes (see Martinez & Saar 2002 



for a description of these processes). It will be shown that 
a suitable shift of the points will reduce 7. However, the 
shift may also be too big and completely destroy the sin- 
gularity of £(r). This is the case, for example, if the shift 
is Gaussian distributed. 

In case of the Voronoi vertices process only simulation 
based investigations are possible. A segment Cox process 
is also studied by this means. The analysis performed on 
both point processes suggest that the effect on 7 is similar 
to that for Neyman-Scott processes. 

2. Point field models 

This section presents some point field models whose pair 
correlation functions have a singularity at r — 0. All of 
them have considerably deficiencies with regard either to 
the exponent in the power-law (Q) or to the underlying 

physics. 

Buryak & Doroshkevich (1996) suggested a cluster pro- 
cess, where the galaxies are situated on randomly oriented 
lines and planes of infinite extent. It models the observed 
fact that galaxies often seem to be situated on filamen- 
tary and planar structures (Bhavsar & Ling 1988). The 
cluster points form one-dimensional Poisson processes on 
the lines and two-dimensional Poisson processes on the 
planes. The two-point correlation function of this point 
field is proportional to r~ 2 for small r, thus 7 = 2. 



A similar process has been introduced by Pons- 



Borderia et al. (1999 ) for testing the performance of dif- 
ferent estimators of the two-point correlation function. A 
segment Cox point field was used for this purpose. The 
process is generated by first placing randomly, with a given 
density X s , segments of fixed length £, within the window. 
Secondly, points are scattered randomly on the segments 
with a given value of the mean number of points per unit 
length on a segment, A^. Stoyan et al. (1995) have cal- 



culated the two-point correlation function of this model 
which reads 



£(r) = 



1 



1 



2irr 2 XJ 2irr\J 2 



(2) 
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Again, a value of 7 = 2 describes the point process 
built by the vertices of a Poisson-Voronoi tessellation, see 
Icke & Van de Weygaert (1987) and Van dc Wcygaert & 



Icke (1 389). The assumption that 7 is equal to 2 has been 
verified by Heinrich et al. (1998| ). 

In all these models 7 is too large. A rather simple point 
field model which may have any 7 between and 3 is pre- 



sented by Stoyan (1994). The points of a homogeneous 
Poisson process of intensity X p define the centres of inde- 
pendent clusters. Each cluster consists of two points. The 
two-point correlation function is 



ft) 



8nX p r 2 ' 



(3) 



where f(r) is the probability density function of the dis- 
tance between the two cluster points. This point field is 
a me mbe r of the class of Neyman-Scott processes, see 
Sect. |3.l[ The size of 7 depends on the asymptotic be- 
haviour of f(r) for r — > 0. If f(r) behaves as r a , 7 is equal 
to 2 — a. Therefore, a suitable choice of f(r) leads to any 
desired order. Unfortunately, this point field is very primi- 
tive and has no physical motivation. Thus, it is unsuitable 
to model the galaxy point pattern and is only of techni- 
cal value. An interesting generalization of this model -the 
so-called Ga uss-Poisson pro c ess- has been presente d and 
analyzed by Kerscher (2001 ). Kerscher et al. (2001 ) have 
shown that this model cannot describe well the distribu- 
tion of the galaxy clusters. 

All of the point fields investigated until now have con- 
siderable deficiencies. Those which have the right 7 are 
quite unrealistic models for the galaxy point pattern. And 
those which are appropriate in some astronomical sense 
have the wrong 7. The next section provides a method 
which yields physically motivated models with a two-point 
correlation function which behaves as desired. 



3. The effect of independent random shifts 

3.1. Arbitrary shifts applied to Neyman-Scott 
processes 

This subsection deals with a class of relatively simple 
Poisson c luster processes, the s o-called Neyman-Scott 
processes Neyman fc Scott (1958). They result from ho- 



mogeneous clustering applied to a Poisson process. The 
points of a homogeneous Poisson process of intensity A p , 
called parent points, define the centres of clusters. The 
clusters are independently and identically distributed. A 
cluster consists of a random number of independently and 
identically distributed daughter points. The random dis- 
tance between a daughter point and its corresponding par- 
ent point is distributed according to a probability den- 
sity function d(r). The two-point correlation function of 
such a process is proportional to f(r)/r 2 , see, for example, 
Stoyan et al. (1995). Here, f(r) is the probability density 
function of the distance between two randomly selected 
daughter points of the same cluster. 




Fig. 1. The quantity F*(r), demonstrated for the two- 
dimensional case. The initial distance between the two 
cluster points c\ and C2 is x, the length of the shift is z. 
The distance between the shifted and unshifted points is 
smaller than or equal to r when the shifted point belongs 
to the dashed line. 

Given such a Neyman-Scott process, we now apply in- 
dependent isotropic random shifts to all its points. The 
isotropic random shift consists of shifting a point by a 
random length in a uniformly distributed direction. The 
random length of the shift has a probability density func- 
tion d*(r). For example, in case of Gaussian shifts d*(r) is 
the Weibull density function with p = 2 and b = l/(2er 2 ) 



• exp 



2ct : 



r > 



(4) 



(this is the probability density function of the length of a 
Gaussian distributed vector in 3-D). 

The random shift changes the distribution of the in- 
terpoint distance between two points of the same clus- 
ter. We denote the corresponding distribution function by 
F* (r) , and the probability density function by /* (r) . The 
resulting two-point correlation function is proportional to 
T(r)/r 2 . 

The aim of the following is to determine the asymptotic 
behaviour of f*(r) for r — -> 0. This will yield a recipe how 
to choose the function d* (r) in order to obtain any desired 

7- 

For a first step, assume that every cluster contains ex- 
actly two points and that only one of them is shifted. 
F*(r) is then the probability that the distance between 
both points is smaller than or equal to r. For example, let 
the initial distance between both points be x and let the 
shift of length be z, see Fig. Then, F*(r) is equal to 
the fraction of that part of a sphere of radius z centred on 
the shifted point, which is not farther from the unshifted 
point than r. 

As the initial distance between the two cluster points is 
distributed according to the probability density function 
f(r) and as the length of the shift is distributed by d*(r), 
we can write 



F*(r) 



md*(z)^^d X dz, 



(5) 



4 



Martin Snethlage et al.: Point field models for the galaxy point pattern 



where A(x; z; r) is the two dimensional volume of the in- 
tersection of a ball of radius r with a sphere of radius z, 
whose centres are at the distance x: 



A(x;z;r) 

4ttz 2 , 



x 




( r 2 _ (a . _ z )2) 



if x < r — z and r > z 
if r — z|<x<?' + 2 
otherwise 



(6) 



(the dashed line in Fig. [j]). 

It is rather complicated to find the integral (||) ; we de- 
scribe it in the Appendix. The result is that the behaviour 
of /* (r) for small r depends only on the corresponding be- 
haviours of d*(r) and f{r). If d*(r) behaves as r a and f(r) 
as r 2-7 (where 7 is the order of the singularity of the cor- 
relation function before the shifts), then f*(r) behaves as 
r 3 ~ 7+a , i. e., the order of the singularity of £(r) after the 
shifts is 7 — (a + 1). It has been reduced, as a + 1 > 
(else, d*{r) cannot be a density function). 

Now, let the number of daughter points per cluster be 
arbitrary and let all of them be randomly shifted. As the 
daughter points are independently distributed, this situa- 
tion is principally equal to that one before, only with an 
initial f(r) which behaves as r 3_7+a instead of r 2-7 (that 
is the resulting density function after randomly shifting 
the first points). It follows immediately that f*{r) be- 
haves as r 3 -7+Q+(Q+i) — r 4-7+2o^ xhis corresponds to 
a singularity of £(r) of order of 7 — 2(a + 1). 

Thus, power-law distributed random shifts with den- 
sity d*(r) ~ r a (— 1 < a < 0) will reduce the order of the 
pole of the correlation function by 2(a + 1). 

3.2. Shifting simulated point processes 

The Neyman-Scott processes are too simple to describe 
the observed galaxy distribution. As we have not found 
analytic expression for the probability density functions 
for more realistic point processes, we have studied the ef- 
fect of random shifts on such processes by means of sim- 
ulations. 

The first model to analyze is the Cox segment process 
introduced in Sect. |[ The realization of the model has 
been simulated in a box of side L = 100 with parameters 
A s = IO- 3 , X e = 0.6, and I = 10 (see Fig. §). 

The second model i s based on Voronoi tessellations . 
Van de Weygaert (199l|); |Van de Weygaert fc Icke (1989| ); 



Icke fc Van de Weygaert (1987 ) suggested that clusters 
of galaxies, the peaks of the matter density field, are situ- 
ated at the vertices of the tessellation. It is remarkable how 
well this model describes the cluster distribution. When 
nothing else is known, one can assume a completely ran- 
dom distribution of the initial points, the nuclei that will 
become the centers of the Voronoi cells, i. e., the empty re- 
gions. Since the initial points form a homogeneous Poisson 
process, the tessellation is usually referred as a Poisson- 
Voronoi tessellation. Fig. |^ shows a realization of this pro- 
cess consisting of 10085 vertices resulting from the Voronoi 
tessellation of 1500 Poissonian distributed nuclei within a 
box of side 100. 




Fig. 2. The upper panels show a realization of the segment 
Cox process. The bottom panels show the Hammer-Aitoff 
projection of the points lying within the maximal sphere 
included in the cube. 



The determination of the two-point correlation func- 
tion of the corresponding vertices process is very compli- 
cated, pleinrich et al. (1998| ) have shown that the singu- 
larity of £ (r) at r = is of order 2. As studies of the real 
galaxy point distribution suggest an order between 1.5 and 
1.9, the vertices process cannot strictly be an appropriate 
model for the galaxy distribution. 

A quite natural way to compensate for this difference 
between the model and reality is to add perturbations to 
each point of the Voronoi galaxy model, as we did for the 
Neyman-Scott models in the previous section. That means 
that, after realizing a sample of this model, each point is 
independently and randomly shifted. This has the addi- 
tional effect that the close attraction between the points 
is relaxed. As it seems to be plausible that the degree of 
the attraction between the points corresponds to the order 
of the singularity, the relaxation should lead to a reduction 
of this order. 

Certainly, there is no prior theoretical distribution for 
the random shifts. But the results of the previous section 
suggest that a power-law distribution of shifts should work 
well. There we determined the effect of arbitrary random 
shifts on the order of the singularity of £(r) in case of 
Neyman-Scott processes. As these processes are based on 
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Fig. 3. The upper panels show a realization of the a 
Poisson-Voronoi vertices process. The bottom panels show 
the Hammer-Aitoff projection of the points lying within 
the maximal sphere included in the cube. 

a homogeneous Poisson process as well as both the Cox 
segment processes and the vertices process of a Poisson- 
Voronoi tessellation, it is plausible to assume that in all 
cases the effect of random shifts is similar. Thus, also in 
case of the simulated processes, any suitable density func- 
tion d*(r) should lead to a reduction of the order of the 
singularity. 

We study this conjecture by simulation. We started 
with the simulated Cox segment process and the Poisson- 
Voronoi vertices process and applied random shifts. The 
random shifts used are distributed according to the prob- 
ability density function 



d*(r) 



0.25r 




-0.75 



for < r < 1, 
otherwise. 



(7) 



The effects of the power-law shifts are demonstrated in 
Fig. U (right panels). 

For comparison, we also tried shifting the points by 
Gaussian-distribution shifts. As Gaussian distribution is 
commonly used when modelling an unknown random ef- 
fect, it can serve as a useful reference model. Thus, given 
the simulated point fields, both the Cox segment pro- 
cess and the Poisson-Voronoi vertices process, every point 
is independently shifted by a three-dimensional Gaussian 



$ « 
♦ ° 



♦ 



♦a? 
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dU 



Fig. 4. Shifting Cox and Voronoi processes. The upper 
panels show the orthogonal projection on the bottom of a 
small parallelepiped, with dimensions 15 x 15 x 20, drawn 
from the Cox segment process (Fig. ||), the lower panels 
show the projection on the bottom of a parallelepiped with 
dimensions 15 x 15 x 10, drawn from the Voronoi vertices 
process (Fig. ||). The left panels demonstrate Gaussian 
shifting, the right panels demonstrate power-law shifting. 
The solid symbols represent the original process, the open 
symbols show the shifted points. 

distributed vector with a — 0.5. Fig. |] (left panels) shows 
the result of this procedure. 

Correlation functions amplitudes usually demand av- 
eraging over a rather large number of pair distances and 
are not easy to use near the pole. The expected power-law 
behaviour can be analyzed better using the correlation 
integral - the so-called if-function, 



K(r) = [ 4tts 2 [1 +£(s)]ds. 
Jo 



(8) 



If K(r) behaves as a power law, we can define the 
correlation dimension D 2 as the exponent of the relation 
K(r) oc r° 2 . 

We used Ripley's estimator for a cuboidal window 
([Baddeley et al. 1993Q: 



K{r) 



V 
~N~ 2 



N N 

EE 

i=l i=i 



Xil) 



(9) 



where V is the volume of the cube, N is the total num- 
ber of points, #(•) is the Heaviside's step function and the 
weights u>ij account for the edge correction, being the con- 
ditional probability that point j would be observed given 
that it lies at distance r from point i. 



6 



Martin Snethlage et al.: Point field models for the galaxy point pattern 



If we find all the distances = |x» — Xj| and the 
corresponding weights Wy- , we can order them by increas- 
ing distance, and define the ordered set {r a }^ =1 with the 
attached weights {w a }a=ii then the estimator ([)]) can be 
rewritten as 



V 



K(r a < r < r a +i) = -rpr 



1 



iV 2 u(r p ) ' 



(10) 



for even indices a, since each r a value is repeated twice. 
This is the sample estimate of the cumulative probability 
distribution function of pairwise separations. For a large 
point sample the function ([l0]) demands a huge memory 
storage. In practice we can approximate this function by 
logarithmic discretization of the total pairwise distance 
interval. 

The results are shown in Fig. |5|. The JT-function is 
shown as a staircase, and the straight lines correspond 
to the expected correlation dimensions (£>2 = 1 for pure 
samples, Z?2 = 3 — 7 + 2(1 + a) = 1.5 for shifted samples). 
We can clearly see that D2 is extremely close to the ex- 
pected value. The small differences between the expected 
power law and the real K{r) dependence are probably due 
to the finiteness of the simulated samples and the algo- 
rithmic (precision) problems in modelling singular point 
distributions. The correlation integral for the Gaussian- 
shifted samples does not reach near the pole and has the 
(Poisson) correlation dimension D2 ~ 3, indicating that 
the correlation structure has been destroyed. 

We have also calculated directly the pair correlation 
functions: these do not reach to the smallest pair dis- 
tances, as explained above. Fig. |^ shows the pair corre- 
lation functions g(r) = £(r) + 1 for the Cox process and 
the Voronoi vertices process, both before and after the ran- 
dom shifts. We have tried different estimators (see Pons 



Borden'a et al. 1999), but no significant differences were 
found amongst them. The correlation functions given in 
Fig. p| were calculated using the Rivolo estimator (Rivolo 
1986y~]As in the previous figure, the straight lines corre- 
spond to the expected order of singularity 7 = 1.5. We 
can see that 7 lies in the expected range, supporting the 
conjecture that the reduction of the order of the singular- 
ity by 2(1 + a), using random shifts with the power-law 
density d* (r) ~ d a is a universal phenomenon. 

We also see that Gaussian shifts have completely de- 
stroyed the singularity. 

4. Summary and conclusions 

1. This paper is based on the empirical evidence that the 
two-point correlation function £(r) of the galaxy point 
field follows for small r a power-law £(r) oc r -7 with 
1.5 < 7 < 1.9. Nevertheless, physically motivated pub- 
lished models for the galaxy point pattern do not have 
this property. 

2. A reason for this dilemma may be that these models 
are usually based on some simplifying assumptions. As 
these ideal conditions are only an approximation of 



reality, it is not surprising that also the point field 
model is an approximation. 

3. One way to cover the difference between theory and 
reality is to add a further random element to the the- 
oretical point field model, i. e. to randomly shift every 
point. 

4. The determination of the effect of randomly shifting on 
the order of the singularity of £(r) is too complicated 
for the general case. Therefore, this paper has inves- 
tigated this question only for a relatively simple class 
of Poisson cluster processes, the Neyman-Scott pro- 
cesses. In this case, there is an easy rule how to choose 
the distribution of the random shifts to get any desired 
reduction of the order of the singularity. If d*(r), the 
distribution of the length of the isotropic shifts, be- 
haves for small r as r a , then the order is reduced by 
2(1 + a). Unfortunately, Neyman-Scott processes are 
not variable enough to model the galaxy point pattern. 

5. The vertices process of a Poisson- Voronoi tessellation 
and the Cox segment process represent the galaxy clus- 
tering better, but their two-point correlation functions 
have singularities of order 7 = 2. This is too large. 
Simulation studies support the conjecture that ran- 
domly shifting the vertices by power-law distributed 
random shifts leads to a reduction of 7 in the same 
way as it does in the case of Neyman-Scott processes. 

6. Thus, a suitably shifted vertices process of a Poisson- 
Voronoi tessellation or of the Cox segment process may 
be an appropriate model for the galaxy point pattern, 
as it is physically convincing and agrees with observa- 
tions. 

7. The method of randomly shifting may be applicable 
also in case of other initial point fields, apart from the 
processes studied here. This may lead to a number of 
reasonable models for the galaxy point pattern. 
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Fig. 5. Effect of shifting on the correlation integral K(r). Label U stands for the unshifted original processes, label 
PL_S stands for the power-law shifted processes, and label G_S stands for the Gaussian shifted processes. The left 
panel shows the Cox family of processes, the right panel shows the Voronoi family. The straight lines represent the 
expected power-law behaviour. 
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Fig. 6. Effect of shifting on the correlation function g(r). The left panel shows the Cox family of processes, the 
right panel shows the Voronoi family. The correlation functions for the original processes are shown by solid bullets, 
the power-law shifted processes by open diamonds, and for Gaussian-shifted processes by crosses. The straight lines 
represent the expected power-law behaviour. 
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Appendix A: Derivation of Eq. (|5]) 

Here, the asymptotic behaviour of /*(r) for r — > is de- 
termined. The starting point is the integral (m 



F*(r) 



f(x)d*(z) A ( X ' ,Z i T ' dxd:. 



4irz 2 



(A.l) 



where /(r) behaves as r 2 ~ 7 (i. e. the order of the singular- 
ity of g(r) before the random shifts is 7) and d*(r) behaves 
as r Q , 7 < 3 and a > —1. 

Elementary considerations show that the asymptotic 
behaviour of F*(r) is equivalent to the case where f(r) = 
r 2 -^l [M] (r) and d*(r) = r a l [0<1] (r), where l [aj6] (r) is the 
indicator function of the interval [a, 6], 



1 for r £ [a, b] , 
for r ^ [a, b] . 



(A.2) 



Thus, the calculations are carried out only for this case. 
We can write 



r r-\-Z 



F*(r) = 



-(x- zf 

l-l-Q 



Axi^z 



dxdz 



T—» 
1 z+r 



-(x- zf 
4xT- 1 z 1 - Q 



■dxdz 



r z — r 
r r—z 



u 2 -~<z a dxdz, 



(A.3) 





and its derivative 

r r+z 

f*(r) = 



2x7-1,2 



1 -,1-a 



dxdz + 



r 3 — f+a 



r(l+a)r(3- 7 ) ,- 7+a 
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2a;T- 1 z 1 
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(A.4) 



where 2F1 is the Gaussian hypergeometrical function (see 



Kratzer & Franz 1960). 



For small r it behaves as 
2 Fi(7-2,a,l + a,-l/r) 

r(1 + ?(7 ( -2) 2 ' Q) ra ' if 2-7 + «<0, 
a r 7-2 



2-7 + a' 



if 2 - 7 + a > 0, 



(A.5) 



and 



2 Fi(3 -7,1-0,4-7,1-1/0 

r(4-7)r(7-2-q) 3-7 - f o _l / n 
— I ±— — ' ) 3 ' r a 7 , if 2 - 7 + a < 0, 



r(l-c^(l-r) 



(2-7 + a)(l-r) 



if 2 - 7 + a > 0. 



(A.6) 



When 2 — 7 + a < 0, all terms behave as r 3 ^7+<\ When 
2 — 7 + a > 0, all linear terms cancel each other out. 
Thus, in any case there remain only terms which behave 
as r 3_7 + Q an d so does f*(r). 
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